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Bone filled with marrow is modeled as poroelastic bore filled with a fluid in the frame 
work of Biot’s theory in wave propagation phenomena. Axially symmetric unattentu- 
ated waves are considered. Frequency equations are derived in the cases of permeable 
boundary and impermeable boundary, and are found to be dispersive. The existing data 
of bony materials has been exploited. Phase velocity versus ratio of wavelength to bone 
diameter curves are plotted. From this model and analytical results, some conclusions 
are drawn. 
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1. Introduction 


Osteoarthritis (OA) affects more than half of the population above the age of 65 [1, 
2] and has significant negative impact on the individuals quality of life [3]. Although 
OA has been considered to be primarily articular cartilage (interface between two 
bones) disorder, this disease is accompanied by change in trabecular bone (compact 
part of the bone) and marrow [1, 2]. Interactions between the bone and marrow 
may be useful in detecting the OA progression. Since direct investigation of bone 
is a destructive one, a study of wave propagation in bony materials is feasible, an 
alternative non-destructive way. It is clear that most of current knowledge of bone 
is based on waves in it. 

Adult human skeleton consists of 206 distinct bones. Each consists of two parts: 
one is cancellous (spongy type, more porous) and other is compact as far as their 
mechanical behavior is concerned. They differ in mass densities. In either case, 
bones tend to exhibit a linearly elastic behavior over small strains, Hooke’s law thus 
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applies. It was suggested that the theory of two phase poroelastic materials rather 
than elastic materials is to be applied in the sense that osseous tissue is considered 
as perfectly elastic and fluid substance filling the pores as a compressible fluid [4]. 
In general, bony materials are non-homogeneous, anisotropic and time dependent, 
however, from the Mathematics point of view, they are assumed as homogeneous 
isotropic materials so that Biot’s theory of poroelastic materials can be used [4]. 

From the Mathematics perspective of poroelastic materials, plane strain vibra- 
tions of thick-walled hollow poroelastic cylinder is studied [5] following Biot’s theory 
of wave propagation in poroelastic materials [6]. Wave propagation in cylindrical 
empty poroelastic bore is investigated [7]. Also, cylindrical stress waves in poroe- 
lastic flat slabs has been studied [8]. With regard to the actual problem considered 
here, we facilitate the Mathematical side of the problem by modelling the bone 
filled with marrow as a cylindrical poroelastic bore of infinite extent containing a 
fluid in the frame work of Biot’s theory [6]. Waves of axial symmetry and unatten- 
tuated are considered here following the paper [9]. We are dealing essentially with 
the interaction between marrow and bone in wave propagation phenomena along 
the axis of bone. Nowinski and Davis studied longitudinal waves in cylindrical bone 
elements and computed some basic parameters of bone [4]. But in this paper mar- 
row is not taken into consideration. Hence the study of bone filled with marrow is 
warranted. Because of axial symmetry, there are two stress components Orr, 7; and 
fluid pressure s. Three boundary conditions on the curved surface pertaining to the 
absence of external loads and permeability of the surface give frequency equations. 
Phase velocity is computed as a function of ratio of wavelength to diameter of bone 
both for permeable boundary and impermeable boundary. 

The rest of the paper is organized as follows. Empty cylindrical bone is consid- 
ered in Sec. 2. In Sec. 3, first behavior of waves in marrow fluid independently from 
the solid bone is investigated; next the bone containing marrow is examined. The 
non-dimensionalisation as well as numerical results are discussed in Sec. 4. Finally, 
concluding remarks are given in Sec. 5. 


2. Empty Cylindrical Bone 


Consider cylindrical polar coordinate system (r,0,z) with z-axis along the axis of 
cylindrical poroelastic bone of infinite extent. Let the bone be homogeneous and 
isotropic. The equations of motion of a poroelastic solid given by Biot [6] in presence 
of dissipation (b) which in terms of displacement vectors are: 


o? 3) 
NV?U +VI(A+ N)je+ Qd = (pune + paj bu — U), 
ə? — F 03 > 
V[Qe + Re = pa (pz u + p22U) — ba Cu — U). (2.1) 


In (2.1), P (= A + 2N), N, Q, R are all poroelastic constants, V? is the Laplacian 


? 


operator, pij are mass coefficients, e and e are solid dilatation and fluid dilatation, 
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respectively. The solid displacement components (u,0,w) which are functions of r, 
z and t (here the problem is plane strain and is independent of 0) which can readily 
be evaluated from the field equations of Biot [6] in matrix notation are 


u 
Cy 

—cos(kz) _ pKı(pr) qK2(qr) kK acar) eit (2.2) 

a kKo(pr) kKo(gr) dKo(dr) l 

— sin( kz) 


where w is the frequency of wave, k is the wavenumber, C1, C2, A; are all constants, 
K,,(«) is the modified Bessel function of second kind of order n, and 


p=k1-@)?, q=h1-G)*, d=kO-)#, G= p> =1,23). 
(2.3) 


In Eq. (2.3), Vi (i = 1,2,3) are dilatational wave velocities of first and second 
kind and shear wave velocity, respectively. Using these displacements into stress- 
displacement relations [6], the relevant stresses (o;;) and liquid pressure (s) are 
obtained as follows: 


Crp Aj, cos(kz) Aig cos(kz) A13 cos(kz) Ci 
Orz | = | A sin(kz) Aəsin(kz) Aə3sin(kz) Co | t, (2.4) 
S A31 cos(kz) Azo cos(kz) 0 Aj 


where 


An = 22K (pr) +{2Np? + (P—2N +67) (p? —k?) + (Q+.R5?)(p? —k?)} Ko(pr), 
As = “Ki (dr) +2NkdKo(dr), 
Ao, = 2NkpKi(pr), 
A23 = N(d? — k?)Ky(dr), 
A31 = (Q + Rô?) (p° — k*) Ko(pr), 
A33 = 0, 
2 —2 

j= E a g 
mi = pu — ibp™t, mi2 = pig + ibp™ t, mag = p22 — ibp™' 
A12, A22, A32 are similar expressions as A11, A21, A31 with p and 6? replaced by q 
and 63, respectively. 


(2.5) 


The boundary conditions to be satisfied on the curved surface r = a to be stress 
free are 


Orr +8=0, orz =0, s=0 (for permeable boundary), 


and 
Orr +8=0, Or =0, 2s =0 (for impermeable boundary). (2.6) 
r 


Equations (2.4) and (2.6) together give a system of three homogeneous equations for 
the constants C1, C2 and A, each for a permeable boundary and an impermeable 
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boundary. In order to obtain a nontrivial solution of this system the coefficients 
matrix must be singular. Accordingly, the frequency equation in the case of perme- 
able boundary is given by 


ayl =0, i, j= 1,2,3, 2.7 
y] 


where aj; are same as A, in (2.4) but are to be evaluated at r = a. 
In the case of an impermeable boundary, the frequency equation is given by 


[byl =0, i,j = 1,2,3, (2.8) 
where 
by = ay 4=1,2 7. = 1,2; 3, 


bs1 = (Q + R6?)p(p? — k?) 1 (pa), 
b32 = (Q ic R63)q(q? —= k?) Ki (qa) and b33 = 0. 


These frequency Eqs. (2.7) and (2.8) will be examined numerically in Sec. 4. 


3. Bone Filled With Marrow 


In this section, first, the behavior of waves in marrow is studied independently from 
the solid part of bone. For axially symmetric waves, the displacement potential 
function @ satisfies the wave equation 


Od 18$ o 184 


zo t-a tio = SS 3.1 
ðr? r Or Oz? ġa? (3.1) 

where cp is the velocity of the sound in marrow. The solution @ is given by 
$ = Jo[rk(C2 — 1)?] cos(kz)e* for ¢>1, 9) 


$ = Ip[rk(1 — C?)2] cos(kz)e* for Ç <1. 


In (3.2), Jo is the Bessel function of zero order of the art kind and I is the modified 
zero order Bessel function of the first kind, and ¢ = z= is the phase velocity along 
z-direction. The solution pertaining to ¢ > 1 sone to conical waves which 
are reflected at the boundary while other one corresponds to Stonely waves that 
occur at the fluid-solid interface. Thus the fluid pressure Py (= -pr $) and the 
radial displacement D, (= a) can be computed. On the boundary r = a, one can 
have 


Pi ____w*psJola(c? = 1)]? p F 
D K@-FAlkae-pje “*S** 2" 
ii 2 a for <1. (3.3) 


Dr KE- Ihla- P)? 
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If the fluid is contained in a rigid wall, that is D, = 0, then the former case 
gives 


2. Oty? 
at P =], (3.4) 
which is an ellipse. 
The later case gives 
2 Gt? 
e "7 Si (3.5) 


which is hyperbola. 

In Eqs. (3.4) and (3.5), a;, 6i are the roots of Bessel functions of order one of 
first kind and second kind, respectively, and D = 2a, diameter of the bone. 

Now we consider bone filled with marrow. The boundary conditions to be sat- 
isfied on the curved surface r = a to be stress free are 


PF =F 
Cra Lor: =0, s=0 (for permeable boundary), 
u D; 
and 
Tr =F 0 s 
aoe ZË on =0, 30 (for impermeable boundary). (3.6) 
u De Or 


The first boundary condition in each case represents a condition at the fluid-solid 
interface, while the other conditions remain the same as in earlier section. Above 
boundary conditions also can be put in the following form 


Or +s=—Ps, Or =0, s=0,u = D, 


and 


However, former form is parsimonious, in the sense, it gives set of three homoge- 
neous equations in three arbitrary constants C1, C2 and A1, whereas, later case 
gives four equations in four arbitrary constants. For a nontrivial solutions of these 
constants the matrix of coefficients must be singular. Accordingly, the frequency 
equations in the case of permeable boundary is: 


ley] =0, ij=1,2,3 for ¢>1, (3.7) 
where 
cu = Ayik(C? — 1)2.4,(ka(¢? — 1)2) + pK (pa) ppw?.Jo(ka(¢? — 1)*), 
cig = Ayok(C? — 1)2 Jı (ka(C? — 1)2) + ¢K1(qa)p pw? Jo(ka(C? — 1)2), 
cis = Arsk(G? — 1)? (ha(¢? — 1)2) + kK (da) pp? Jo(ka(¢? — 1)2), 
Cpe Ag t=; J = 1,25. 
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For the case Å < 1, the frequency equation is given by 
ldy|=0, i, j= 1,2,3, (3.8) 


where dj; are same as cij with (C? — 1)? is replaced by (1 — ¢?)?2. 
Similarly, the frequency equations can be obtained for the case of impermeable 
boundary. All these frequency equations will be examined numerically in Sec. 4. 


4. Numerical Results 


In absence of dissipation, the frequency equations obtained in last section are anal- 
ysed by introducing the following non-dimensional parameters: 


P Q R N 
a=—, @=—, a=—, a=— 
1H? °F H ATH Oe 
Fie Pi, i= ay = p22 
p p p 
i Vi? V? W\? 
ee ees =; ce = [| =). 4.1 
+= (E). (E) E). (4.1 
2 H 
p = pu + 2p12 + p22, H = P++2Q +R, r 


>|2 


C : : : Ww 
m = — non-dimensional phase velocity, c= T and cå 


Using these non-dimensional parameters into the frequency Eqs. (2.7) and (2.8) of 
empty bone, one obtains an implicit relation between phase velocity (m) and the 
ratio of wavelength to diameter (4). Phase velocity is computed against this ratio 
both for permeable boundary and impermeable boundary. These equations yield 
infinitely many values for phase velocity, of which, first three values are presented 
in Fig. 1. The values of bone poroelastic parameters A, N, Q, R and its mass coef- 
ficients p;; are computed as in the paper [4]. The values of Young’s modulus and 
Poisson ratio are taken to be 3 x 10° and 0.28, respectively as suggested in the paper 
[4]. Velocities of dilatational waves and shear wave (Vi, V2 and V3) are computed 
using Biot’s theory [6]. From Fig. 1, it is clear that the values of permeable boundary 
and the values of impermeable boundary are almost equal except in the neighbor- 
hood of 2 = 1.5. In this neighborhood, the values of impermeable boundary are 
greater than that of permeable boundary. This deviation decreases from the value 
1 to value 3. The first values of phase velocity for permeable boundary are almost 
same. In order to investigate the effectiveness of poroelastic constants A, N, Q, R 
on phase velocity, we make sensitivity analysis by changing the Poisson ratio. It 
is clear from the paper [4] that poroelastic constants A, N,Q, R are functions of 
Poisson ratio. Therefore, by changing the Poisson ratio values from 0.26 to 0.3 in 
step of 0.1, the pertaining values of A, N,Q, R, Vi, V2 and V3 are computed. Then 
for each set of these values, the phase velocity is computed against (4) both for 
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Phase Velocity (m) 


| 
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Fig. 1. Variation of phase velocity with wavelength/diameter in the case of empty bone for both 
permeable boundary and impermeable boundary. 
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Fig. 2. Variation of phase velocity with wavelength/diameter at various poisson ratio in the case 
of empty bone for permeable boundary. 


permeable boundary and impermeable boundary. The results of permeable bound- 
ary and impermeable boundary are depicted in Figs. 2 and 3, respectively. From 
these figures it is clear that Poisson ratio influences phase velocity. In general, the 
values of permeable boundary are greater than that of impermeable boundary. The 
difference between the values of permeable boundary and impermeable boundary 
increases as Poisson ratio increases. In the case of empty bone, phase velocity is 
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Fig. 3. Variation of phase velocity with wavelength/diameter at various poisson ratio in the case 
of empty bone for an impermeable boundary. 


sensitive with respect to nature of the boundary and Poisson ratio. Using non- 
dimensional parameters into the frequency Eqs. (3.7) and (3.8), one obtains an 
implicit relation between phase velocity (Ç) and the ratio 4 for fixed E and ma 
As in the case of empty bone, sensitivity analysis of phase velocity against Poisson 
ratio is made in the case of bone with marrow. For the numerical work, the veloc- 
ity of sound in marrow (cy) and the density of marrow py are taken to be 93 m/s 
and 900 kg/m3, respectively. These values are obtained in a private communication 
[10]. It is obtained that the phase velocity is not depending on Poisson ratio unlike 
the case of empty bone. Moreover, from the calculations it is seen that values for 
permeable boundary and impermeable boundary are equal. This means that nature 


Phase Velocity (m) 


Ratio of Wavelength to Diameter 


Fig. 4. Variation of phase velocity with wavelength/diameter in the case of bone with marrow. 
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of boundary does not have any influence over the values in the case of bone filled 
with marrow unlike empty bone. All these common results are depicted in Fig. 4. 
From Fig. 4, one can infer that first the values increase as the ratio increases up to 
3 = 1.25, next decreases up to 1.5 and then increases linearly after 1.5. From the 
figures, it is clear the values pertaining to empty bone are greater than that of bone 
filled with marrow. From the figures, it is evident that dispersive behaviors in the 
empty bone and the bone filled with marrow are different in phenomenon, in the 
sense, phase velocity is sensitive with respect to nature of boundary and Poisson 
ratio for empty bone, while it is not in the case of bone filled with marrow. 


5. Conclusion 


Cylindrical Bone filled with marrow is modelled as a cylindrical bore filled with 
fluid in the frame work of Biot’s poroelastic theory to investigate the interaction 
at the interface of bone and marrow. A comparative dispersive study has been 
made between the empty bone and bone filled with marrow for both permeable 
and impermeable boundaries. The available data of bony materials has been used 
for the numerical evaluation. One limitation here is existing data is based on certain 
approximations. For more reliable data of bony materials, one has to make extensive 
experimental work. However, this Mathematical model and numerical evaluation 
give some rough estimate with regard to bone, in turn, this rough estimate may 
provide some deep insight into the mechanics of bone. 
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